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0> ' Abstract 



We measure the diffusion constant for Chern-Simons number for classical, lattice SU(3) 
Yang-Mills theory, using a generalization of the topological definition of Chern-Simons num- 
^ I ber developed recently by Moore and Turok. The diffusion constant is much larger than that 
for SU(2), even before the ratio of coupling constants has been accounted for, which implies 
that chiral quark number is efficiently destroyed by strong processes during the electroweak 
phase transition. For the physical value of ag we estimate the decay time for chiral quark 
^ ! number to be about 80/T, although various systematics make this number uncertain by 
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5 '. 1 Introduction 

Baryon number is violated in the standard model and there has been a growing interest 
in trying to understand mechanisms which might use this violation to generate the baryon 
^ , asymmetry of the universe during the cosmological electroweak phase transition. 

All the needed ingredients for generating a baryon asymmetry [|I|] are present; the violation 
^ ; of baryon number shuts off abruptly while the plasma is out of equilibrium due to the 



about a factor of 2. 

PACS numbers: ll.lO.Wx, 11.15.Ha, ll.15.Kc 



^ ■ motion of a bubble wall (phase boundary), and if there is C and CP violation then these 
conditions can give rise to a net baryon number generation. Interest in this scenario has 
focused especially on a particularly efficient mechanism in which the CP violation, in the 
form of spatially varying Higgs condensate phases on the bubble wall surface, generates a 
chiral top quark asymmetry 0. This can be transported by particle diffusion ^ into 
the symmetric electroweak phase; the left handed quark number then biases SU(2) winding 
number changing transitions ( "sphalerons" ) which generates a net baryon number. 

A complication to this scenario is that chiral quark number is damped by SU(3) color 
winding number changing events |^ , similar to the phenomenon responsible for the spon- 
taneous breaking of chiral symmetry in QCD. To be more concrete, if there is a chemical 
potential fi for chiral quark number, then the total chiral quark number density (left handed 
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quarks minus their antiparticles minus righthanded quarks, plus their antiparticles) is 



12 12 

species 

where the 8 comes from summing over up and down types, particles and antiparticles, and 
the two chiralities, and the T^/12 is the (leading order in fi and g"^) change in number density 
due to the chemical potential for an ultrarelativistic fermionic species. {Nc = 3 is the number 
of colors and Np = 3 is the number of fermion generations.) The free energy liberated by a 
strong sphaleron event is fi for each left handed particle destroyed or righthanded particle 
created, which equals 4:Nf, since in each generation there is a creation of a righthanded 
particle and a destruction of a lefthanded particle, both for up and down type quarks. 
The rate per unit volume of strong sphaleron transitions will then be ANp[ix/T)T ss-, where 
Tasi the linear response coefficient of strong sphalerons to a chemical potential, equals half 
the diffusion constant per unit volume of SU(3) Chern-Simons number, by a fluctuation 
dissipation relation [0, 0. Since each transition changes chiral quark number by —4:Nf, 
one findfl 

^ = (4/V,)^iir„ = -Hl^Q.I^l. (2) 

dt ^ ^> T Nc T3 ^ ' 

The time constant for the decay of is therefore 

Even without calculation we know that is much larger than the corresponding weak 
sphaleron rate Ty^s, because SU(3) contains SU(2) and the strong coupling is larger (so 
nonperturbative physics sets in on a shorter length scale). Hence chiral quark number in the 
symmetric electroweak phase decays mainly through strong phenomena. If the time constant 
is shorter than the typical time a quark reflected from or moving off of the bubble wall spends 
in the symmetric phase before the wall catches up with it, then the strong sphaleron rate will 
be relevant and will reduce the baryon number generation. In this case we may only need 
to know the ratio of the strong and weak sphaleron rates to determine the baryon number 
generated. If r is shorter than the time it takes a particle to get from the middle of the 
bubble wall to the symmetric phase, then the chiral quark number may be destroyed before 
it reaches an environment where it can generate baryon number, and the strong sphalerons 
will qualitatively reduce the production rate for baryon number. (It is important to treat all 
relevant processes, such as the difference in diffusion constants between right and left handed 
quarks, and the conversion between chiral top quark number and Higgs particle abundance 
(net Higgs particle hypercharge) due to the top quark Yukawa coupling [Q; these might 
change this picture somewhat.) In any case it seems that the investigation of F^s is well 
motivated. 

It is conventional to write F in terms of a dimensionless constant k, 

Tss = KaaatTy2, (4) 



^ The derivation here follows that in Q except that they miss the factor of in the relation between 
Qs and /it, so their rate constant for the decay of Q5 is 3 times too large. 
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where the factor of 2 is because in conventional usage Kss denotes the diffusion constant for 
Ncs, not the response coefficient, k^s may depend nontrivially on couphngs and the particle 
content. One should take as at a renormalization point on order the temperature scale. 
Beyond leading order, k must have a logarithmic renormalization point dependence; we will 
discuss the "best" renormalization point in the next section. 

A first attempt to compute F^s was made in 0, using the classical approximation for the 
infrared dynamics first suggested in |jlO|] and a numerical implementation of Yang- Mills theory 



and of Chern-Simons number developed and used in |]TT], |12[. Unfortunately, this definition 
of Chern-Simons number suffers from lattice artifacts, leading to a spurious ultraviolet signal 
and an incorrect normalization of the response to the real, infrared winding number change. 
The lattice implementation of Yang-Mills theory also requires perturbative corrections to 
the tree level match between lattice and physical length scales It is known in the case 



of SU(2) that both corrections are numerically important. The SU(2) diffusion constant 
without these corrections gives a lattice spacing independent k [|12], which is qualitatively 
different from what is expected on theoretical grounds. Arnold, Son, and Yaffe have argued 
that the interaction of the infrared modes responsible for winding number change with "hard" 
short wavelength excitations is essential, and because of it k should depend on the product 
of the plasma frequency ujpi and the inverse nonperturbative length scale /„p oc l/g'^T as 

K oc (Wpi/np)"^ , (5) 



at least when {oJpilnpY ^ 1 13 • Since ujpi on the lattice depends on the lattice spacing a 
as a"^/^ [|^, this implies that k should be proportional to a; and for the physical quantum 
system ujpi ~ gT, so Kyjs oc aw 

Recently a topological technique for measuring winding number change on the lattice has 
been developed for classical SU(2) gauge theory Using it, and applying the corrected 



matching between lattice and physical length scales, reveals that T^s does depend on lat- 
tice spacing, verifying that previous results were contaminated with lattice artifacts. The 
coarsest lattices used there were insufficient to achieve {oJpilnpY ^ and the dependence 
was (therefore?) weaker than linear with a; but we will assume here that the reasoning of 
Arnold, Son, and Yaffe is correct and that for sufficiently large hard thermal loop effects the 
rate does scale with {ojpilnp)~'^ ■ 

Because F depends on the physics of hard modes, which is definitely not reproduced 
correctly on the lattice |]T^, we do not know how precisely to convert F measured there into 
F in the quantum theory. This problem was addressed recently by Arnold, who argues that a 



reasonably accurate conversion is possible [15|. The small a limit of the ratio of F for SU(2) 



and SU(3) should also be insensitive to this problem, since the distortion of the hard thermal 
loops is common to the two lattice theories. Hence we can compute the ratio F^s/F^is on 
the lattice with smaller systematics. The purpose of this letter is to extend the results of 
and the technique of |T^ to SU(3), and to use them to compute F^^ on the lattice, and 



hence to find Vss/'^ws and to estimate F^ 
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2 Diffusion constant: General discussion 



Most analytic work on F has discussed the broken phase of Yang-Mills Higgs (YMH) the- 
ory. In this case the gauge connection is usually close to some topological vacuum, and 
the diffusion rate is controlled by the free energy of configurations midway between vacua, 
"sphalerons." This way of viewing things can also be useful in the symmetric phase or in 
pure Yang-Mills (YM) theory, where the diffusion rate will also depend on how often the 
system is straddling between vacua. 

Consider classical YM or YMH theory in a fixed finite volume and regulated at a length 
scale much smaller than Xjg^T. (We are thinking of lattice regulation, but we will use 
continuum notation here for convenience.) We can define the closest vacuum to the spatial 
connection A at time t by "cooling" [|T^. That is, consider evolving A under (gauge invariant) 



straight dissipative dynamics, 

dA[x, t, t) _ dH{A{t,T)) 
dr dA{x, t, t) 



A{x,t,0) = A{x,t) , (6) 



with H the Hamiltonian. (The cooling time r is not to be confused with the decay rate 
for chiral quark number discussed in the introduction.) Such cooling was recently used as a 



technique to improve the local operator method for tracking Ncs [T^l- At sufficiently large 
r the connection will settle into a vacuum configuration, and we define this to be the nearest 
vacuum. 

It will sometimes occur that there is a time tgph such that the nearest vacuum at time 
tsph + e has a different winding number than the vacuum at time tsph ~ ^-0 At time tgph 
the dissipative evolution will never get to a vacuum state, and just before and after, it will 
take a very long time. We can define the system to be "in a sphaleron" if the cooling time 
r required to get close to vacuum exceeds some threshold rthresh- Here "close to vacuum" 
can be given a rigorous definition, eg the total action of the remaining cooling path to the 
vacuum is less than (vr/^f)^. 

The probability that the system is "in a sphaleron" depends only on Tthresh and the 
thermodynamics of the spatial connections. Further, if we define 

Etix,t,T) = [Do, D,nx,t,T) (7) 

(with Dq the covariant t derivative), then since the distribution of values for E{t, r = 0) is a 
thermodynamic property, and since dE/dr depends only on the connections and on E, |Jl9[ 



dr ' ' 'dAf{x,t,T)dA''j{x,t,T)' ' ' 

then E{t, r) is also distributed according to thermodynamics. If we define the sphaleron 
narrowly enough that the system typically remains "in a sphaleron" for a time short com- 
pared to the inverse plasma frequency, so E{t, 0) does not change much from the beginning 

■^This statement also has a rigorous gauge invariant meaning, even on the lattice, if the lattice spacing is 
not too coarse: the winding number difference of the vacua is defined as the winding number of the gauge 



transformation carrying one to the other, measured using the algorithm of Woit |20 , or, for SU(N), the 
technique developed in the next section. 
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to the end of the sphaleron event, then the length of a sphaleron event will also depend only 
on Tthresh and on thermodynamics, and the total spacetime density of sphaleron events will 
be a thermodynamic property, depending on the thermodynamics of A and E alone. We 
have not shown here that this spacetime density has a good large volume limit, but we will 
assume this to be the case. 

Since the spacetime density of sphaleron events depends only on thermodynamics, we 
know that the quantum theory value is reproduced by the classical theory in the a ^ 
limit, with corrections due to the thermodynamics of the A fields which are O(a^) if we use 
the dimensional reduction calculation |^ to establish the value of the coupling constant of 
the 3-D theory. There may also be 0{a) corrections in the thermodynamics of the E fields, 
ie in the relation between time scales p2[, which have not been calculated. Further, since 
the only length scale in the thermodynamics of the A and E fields is l/{g'^T), the spacetime 
density of sphaleron events can be written as k,i{Nc) x [aT)^, with k,i{Nc) a pure number. In 
YMH theory, ki also depends on the Higgs potential parameters x and y. It should approach 
the YM theory value in the limit of large positive y, deep in the symmetric phase, and it 
becomes exponentially small at large negative y, deep in the broken phase. 

Each sphaleron event changes the winding number of the nearest vacuum by ±1. If the 
signs of each change were independent, then F would equal Ki{Nc)a'^T'^/2. But the signs will 
in general be correlated, F = K2Ki{Nc)a'^T'^/2, where K2 describes the degree of correlation 
in the signs of sphaleron events, and depends on the dynamics. In particular, Arnold, Son, 
and Yaffe argue that plasma oscillations will make the system go back and forth through 
sphalerons of opposite sign. On short time scales the motion of infrared magnetic fields will 
be oscillatory, and on long time scales it will be overdamped. They conclude that the system 
will go through on order {upi/ g'^T)'^ ~ l/{g'^h) sphaleron events per permanent winding 
number change, so K2 ~ g^h fl^. In the classical theory, the role of h is played by the 
regulator scale, and K2 oc g'^aT /A = Pj^^, at least for large Pl jTsj . 

To determine F correctly, it is necessary to count winding number changes correctly; 
to get the thermodynamics, and hence ki, right; and to get the dynamics, and hence K2, 
right. The first two problems are separate from the third, and we deal with them in the next 
section. Getting the dynamics right is harder, and in our opinion this problem has not been 
solved. However, the work of Arnold [|T^ suggests that lattice results can be converted to 
continuum results with fairly modest systematic error, and we will use his matching here. 



3 From SU(2) to SU(iVe) 

There are no complications in extending the standard Kogut-Susskind implementation of 
3+1 dimensional SU(2) Yang Mills theory [^] to SU(3), and the thermalization algorithm 
for the SU(3) case was developed in |^. What remains is to extend the one loop matching 
of the thermodynamics of lattice and continuum systems, and to extend the topological 
tracking of winding number, from SU(2) to SU(3) (or SU(A'c)). 
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3.1 thermodynamics 

We deal first with the thermodynamics. As we discussed in the last section, it is only the 
thermodynamics of the spatial connections (and of the E fields) which are important; so our 
goal is to make sure that the thermodynamics of spatial connections at finite lattice spacing 
a are as close as possible to the continuum thermodynamics. Here and throughout we will 
use the notation of ||12|, We will not attempt to make this section self-contained; the 
reader is referred to [113 for details on the approach. All we do here is generalize to SU(Ai'c) 



the gauge field part of the SU(2) calculation done there. The details of this section are not 
important in what follows, only the final result, so the uninterested reader can skip to the 
next subsection. 

The thermodynamics of the real time system we are considering are determined by the 
path integral 

Z = JvAAoexp{-(3LHL), (9) 
Hl = Y: (^-iReTrf/,,(x: 

x,i<j ^ ^ ^ 

+ E lmx)Ao{x + .)f//(x) - A,{x)f + E ^Alix) , (10) 

X,i X 

with the bare Debye mass m|, = Here Uij{x) is the elementary plaquette which 

extends from x in the i,j directions. The gauge coupling has been absorbed into Pl which 
(at tree level) equals^ = 4:/{g'^aT). We want to improve the Hamiltonian so the thermo- 
dynamics produced by this partition function match more accurately those of the continuum 
system. The idea is that, since the lattice and continuum theories only differ strongly in the 
ultraviolet, one should compute the influence of ultraviolet modes on the infrared physics 
perturbatively, in the lattice and continuum theories, and find the difference (which is free 
of infrared divergences) . Because the infrared length scale is well separated from the length 
scale where the lattice and continuum theories significantly differ, the difference can be writ- 
ten as an operator product expansion, and only the super-renormalizable terms are needed. 
One compensates for these terms by making shifts in the wave functions and couplings of 
the theory, thereby correcting the lattice theory in the infrared for its ultraviolet differences 
from the continuum theory. Because the theory is super-renormalizable, one loop perturba- 
tive corrections are 0{a), leaving the infrared behavior of the lattice and continuum theories 
matching up to O(a^). 

The exception to this rule is dimension 2 operators, where the one loop contribution is 
0(1/ a) in physical units, which in lattice units is ~ /J^^^, and a full 0(a) correction 
requires a three loop calculation. One dimension 2 operator, the Debye mass, appears here. 
However, the large value of m|, makes the influence of the Aq fleld on the thermodynamics of 
the gauge flelds perturbative; to study their thermodynamics, we will want to integrate out 
the Aq fleld, and the 2 and 3 loop corrections to m|, and 1 loop wave function corrections 
to the Aq fleld will only change the result of that integration at order (31'^^'^ , as we discuss 

''For Nc ^ 2 this notation differs from that usually used in 4 dimensional lattice QCD, which uses 
1 - (l/7V)Re Tr. 
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below. Hence, we will not calculate these corrections here. But we will need to know the 
Debye mass. Since the bare value is zero, m|, equals the counterterm, which is computed at 



one loop in [24 



= (lattice) =^<,*r^ (physical). (11) 

It only remains to find the correction to the plaquette term, which must be multiplicative, 
ie the term 1 " (l/2)Trf/ above is multiplied by a wave function correction Z^. This can 
be absorbed into (or understood as) a shift in i3l- Contributions to Za arise both from 
self-energy and vertex corrections, and are actually easiest to compute in the theory with 
Ns fundamental scalars added. The gauge field is properly normalized if the full effect of a 
gauge line propagating between scalar lines is the same in the lattice and continuum theories. 
The scalar wave function receives a renormalization which, by minimal coupling, changes the 
strength of the scalar-gauge vertex. There are also loop corrections to the vertex and to the 
propagator, illustrated in Figure |l|. What enters the calculation is the difference between 
the loop corrections to the gauge-scalar vertex and the scalar wave function, and the pure 
gauge and Aq contributions to the gauge self-energy. The only differences between the SU(2) 
and SU(iVc) calculations are the group factors; no new diagrams or new momentum integrals 
appear. Examining |jl3|], one finds that the fundamental scalar corrections to the self-energy 
depend on NsTiTaTb = {Ns/2)Sab, and that all but one of the other diagrams have group 
factors proportional to Nc (if one takes the appropriate differences between contributions to 
the scalar- vector vertex and the scalar wave function correction). The diagrams needed and 
their values are tabulated there. 

The one exception is a contribution to the gauge field tadpole diagram arising from a 
term in the gauge field 4-point interaction, which comes entirely from anticommutators of 
Lie algebra generators, and has a group factor of 

1/2 \ 

3 ViV ^^'^''^'"^ ~^ ^ac^bd + ^ad^bc] + dabedcde + dacedbde + dadedbcej (12) 

and a Lorentz dependence of [1/2) F^j- This term is purely nonrenormalizable operator 
and has no analog in the continuum theory; it also dominates the correction between lattice 
and continuum theories. To find the group factor for the contribution to one contracts 
against 6cd] the 6ab type terms give (2/3)(A^+ 1/A^) and the dabe type terms give (2/3)(A^ — 
4:/N). The space integral gives 1/(3/?^). 

All told, the contributions to the gauge field renormalization are (^ = 0.152859, S = 
3.17591) 

PL{ZA-l) = l(2N^-i-]+Nj^-l]+Nj^-f]+Nj^--L] , (13) 



9V NJ \127T 9 J \247T 6n J \48n 12tx ^ 

where the first term is from the contribution discussed above, the second term is from all 
other gauge field contributions, the third is from the (adjoint scalar) Aq field contribution to 
the self-energy, and the last term is if there are fundamental scalars present. If one puts in 
a "naive" bare value of PiMre then the simulation is equivalent to one with the appropriate 
Za and PL,imp = Z^^ jSifiare-, which numerically is 

/AN 2 \ 

(3L,imp = PiMre ' (^-^ - ^ j - .0389iV, - .0340iV, - .0170iV, (14) 
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Figure 1: Examples of diagrams contributing to the renormalization of the gauge field. The 
full effect of a gauge line propagating between scalar lines must match between theories; be- 
sides self-energy corrections, there are also corrections at the vertex, and corrections because 
the vertex is renormalized at the same time the scalar propagator is. 



where again the contributions are from the anticommuting part of the tadpole diagram, other 
vector contributions, the Aq field, and any scalars present. For A^^ = 3 and no fundamental 
scalars, we find l3L,imp = PL,bare — 1.330. This gives us the correction we need to convert from 
lattice to physical length scales. 

Two comments are in order here. First, the correction is totally dominated by the 
tadpole term. This behavior is typical in lattice gauge theory. It is caused by the compact 
nature of the connections and it is the reason a perturbative matching between the lattice 
and continuum theories is so necessary. Second, the correction is larger for SU(3) than for 
SU(2). Since the diffusion constant depends on the fourth power of the conversion between 
length scales, one must ensure that the matching is quite good. If for instance we assume 
that there are unknown O(a^) tadpole type corrections of magnitude equal the square of 
the 0{a) correction, and we ask that these corrections be at the 1% level (which will still 
give 4% systematic errors), we need Pl > 13; we must work on quite fine lattices. Though 
this increases the numerical demands, it is not all bad, since it means that the elementary 
plaquettes will be quite close to the identity and it should be possible to make the connections 
quite smooth. We will need this property in the next subsection. 

Now we have related the lattice theory to the continuum theory with Ao field, and with 
the wrong Debye mass. We need to understand how the value of mjj modifies the gauge 
field thermodynamics, and this is easiest done by integrating out the Aq field. In terms of 
the natural 3-D length scale 1/g^T, the Debye mass is mj^ ~ Pl^q'^T)^, so it is indeed heavy 



enough to integrate out. At one loop, the modification to the gauge coupling is ET 



This correction is formally 0{(3i^^'^). Now recall that mjj receives two loop corrections 
which are 0(/3^^) in lattice units, or 0{g'^T^) in physical units. One can see by plugging 
m|, = {Af]L + B)g^T^ into Eq. (p!5|) that this correction leads to an 0(/3^^''^) correction to g"^. 
Similarly, the 0(/?^^) correction to the Aq wave function renormalization can be absorbed 
by a rescaling of Aq into a shift in m'j^ of order 0(m|)/3^^), also leading to an 0{(3^^^'^) 
correction. We will neglect these 0(/3^^^^) corrections. 
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In addition to the 0{f3i^^'^) correction we have just mentioned, the integration over the 
Aq field gives an 0(/3^^) two loop correction, and induces nonrenormalizable operators which 
affect physics at the nonperturbative length scale by 0(/3^^^^). While the nonrenormalizable 
operators are ignorable. The two loop correction is parametrically as important as the 0(/?^^) 
correction to the lattice-continuum match which we have just calculated. However, we 
believe that it is numerically much smaller. This is because of the "tadpole" character of the 
dominant 1-loop effects we have studied. It can be argued that the expansion parameter for 
the integration over the Aq field is g'^T / A-nmi:) ~ 0.1/3^^''^, whereas the "tadpole" corrections 
were ~ 1 x (31^. Numerically, at (3l = 16, the one loop, 0{[31^^'^) correction to g"^ from 
integrating out the Aq field, Eq. (0), is less than 1%, while the 0{l3j}) correction in the 
lattice-continuum match is ~ 8%. Hence, although the two loop correction from integrating 
out the Aq field is formally of the same parametric order as corrections we include, it should 
be numerically unimportant. Note that the numerical unimportance of this two loop effect 
is also important in the study of thermodynamic properties of SU(2) YMH theory using the 
dimensional reduction program if one integrates out the Aq field in that program. We will 
only make corrections due to Eq. (|14D and Eq. ([TsD in this paper. 

Finally, there are also 0{a) lattice corrections to Eq. (0), which can be understood as 
corrections in the conversion between lattice and continuum time units. These have not been 
calculated, but it was argued in that they are dominated by tadpole effects which are 
the same as those occurring in the lattice-continuum match for the spatial gauge fields. We 
will use the prescription proposed there, leading to another 0(a) but probably modest error. 

The thermodynamics are now under control, and while there will still be uncorrected 
0(a) errors in ki, the large "tadpole" type corrections are taken care of. 

3.2 winding number 



Now we discuss the extension of the winding number tracking technique of ||T6[ to SU(A''c). 

The idea of that paper is to keep track of a notional group valued scalar field 5*, with 
Hamiltonian 

Hs = Y.^-\^^ '^^S\x)U,{x)S{x + i) . (16) 

One then evolves S dissipatively and agressively to minimize Hs- The Chern-Simons number 
of a configuration with S = I everywhere is approximated to be zero, and the total winding 
number change during an evolution is tracked by gauge transforming to the gauge S = I 
whenever that gauge is everywhere smooth, ie there is no neighborhood where the connection 
matricies U are far from the identity. When the winding number of the underlying gauge field 
configuration changes, then S will go through a period where it is not smooth somewhere, 
as it adjusts to describe the new winding number state. When it has returned to being 
everywhere smooth, the gauge transform to 5* = / is a large (but smooth) transformation; 
we find its winding number and use it to increment the cumulative winding number change 
to date. The winding number of a gauge transformation is determined with an algorithm 
which is essentially that of Woit [pO| . 

Two things become more complicated when one goes to SU(A''c). The first is the imple- 
mentation of the dissipative algorithm for S. The basic element of the dissipative algorithm 
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is to minimize Hs with respect to S at one site x p5|. The easiest way in SU(2) to find 
S{x) which minimizes Hs is to sum the parallel transports of nearest neighbors, which will 
be a constant times the desired element of SU(2), and to project the modulus to SU(2). For 
SU(iVc) the sum of several group elements is not generally a multiple of a group element, 
and one must orthogonally project to SU(A^c) by a more complicated algorithm. First, scale 
the matrix so the modulus of its determinant will be close to 1. Call the resulting matrix 
M. To project towards U(iVc), one repeatedly replaces M (3/2)M - (l/2)MMtM; if this 
process converges it gives the closest element of U(A'"c) to the original matrix M . Then one 
makes a U(l) rotation by the angle (— 1/3) argDetM to get to SU(Ai'c). If the projection 
to U(A'"c) fails to converge, for instance because the slave field is varying wildly around the 
point, then one does not update at this point; but in practice this essentially never happens. 
The algorithm to combine these elementary projections into an efficient quench is the same 
as in the SU(2) case |T^ . 



The second complication is determining the winding number of a slave field configuration. 
For SU(2) there is a simple algorithm because the group has the same dimension as the space. 
Interpolating S{x) between lattice sites by a geodesic rule, S becomes a map from to 
SU(2)= 5''^, which are of the same dimension, and the winding number is just the oriented 



sum of times some fixed point in SU(2) is covered ||2^. There is an extension of this idea 
to SU(3) by choosing a cannonical map of SU(3) with a dimension 4 subspace cut out into 

X such that the part carries the relevant topological information [^. The 4 
dimensional excision will not generally be important since SU(3) is 8 dimensional and we 
are mapping from a 3 dimensional space. 

However, because the underlying fields are quite smooth, it turns out there is an eas- 
ier way to determine the winding number. One constructs the vacuum field obtained by 
gauge transforming the naive U = I vacuum by S, and then estimates its Nqs directly by 
integrating eijkfabcKAlA';. Defining 

A^ix + t/2) = ^ReTr - iX'^Uiix) = ^ReTrS^x + i)(-zA")5(a;) (17) 



and 



the integral is 



^ ^ A-{x + t/2)+Anx-z/2) 



Ncs = 7^Y: fatcAUx)Al{x)A^{x) . (19) 



The result is not gauge invariant and will not be an integer. However, if the slave field is 
smooth, then the result will be close to an integer; if the slave field used in a gauge trans- 
formation is always suitably smooth then it will be possible to unambiguously reconstruct 
what integer the above integral was "trying to give us". We have tested this idea with SU(2), 
where it is possible to compare the value of the integral to the (integer) winding number 
determined topologically. This is illustrated in Figure |^; as seen there, the values from Eq. 
(p!9|) always understate the winding number, but by a fairly consistent amount. The values 
obtained for gauge transformations of one winding number do not overlap those which arise 
from another winding number, so one can unambiguously reconstruct the winding number 
of gauge transformations from the values of Eq. ([I9|) . To implement the same idea in SU(3), 



10 
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-2-1012 -2-1012 

Measured winding numl^er Reconstructed SU{3) winding number 

Figure 2: The value of the integral, Eq. (|19D , versus the actual winding number for gauge 
transformations in an SU(2) evolution at /S^ = 8 on a 20^ lattice, at left, shows that the 
integral can be used to unambiguously reconstruct the original winding number. At right, 
Eq. ( |T9| ) for gauge transformations in an SU(3) evolution on a 24^ grid at (3l = 16, plotted 
against the reconstructed winding number. The values of the integral are clustered with 
large breaks which makes the reconstruction unambiguous, even though we have no direct 
integer measure of the winding number. 

one just writes the time and value of Eq. (|19D every time a gauge transformation is made. 
Plotting the values of the integral which occurred, one finds the appropriate breaks and can 
then reconstruct the winding number changes, also illustrated in Figure ^ There is never 
any difficulty in the reconstruction for the lattice spacing and volumes used in this paper, 
though presumably this technique should start to have problems on very coarse lattices or 
large volumes. 



4 Numerical results 



It appears from theoretical arguments |T^, |TB| and numerical results that the Ncs 
diffusion constant in SU(2), k^s, depends on the lattice spacing, because finer lattices have 
more hard modes contributing to hard thermal loops jl^. Since it is more numerically 
expensive to study SU(3) it makes the most sense to try to come to grips with this problem, 
and with the problem of trying to include hard thermal loops properly, in the SU(2) theory. 
Hence we will determine the lattice value of Tss at only one lattice spacing, albeit a fairly 
fine one. It should be possible to use this value to establish the ratio Tss/T^s, up to errors 
from how F approaches the large ujpi scaling regime. The absolute value of F^^ can also be 
estimated using the arguments in |15|, but the systematic errors will be bigger here. 

It is not known what finite volume systematics may occur in the calculation of F^s, so 
we measure it on a range of (cubic toroidal) lattices. Our results are for PL,bare = 16 and for 
8^, 12^, 16^, 24^, and 32^ lattices. The results are tabulated in Table |^, which presents Kss, 
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4.0 ± 0.8 


16 


16 


8.3 ±1.1 


16 


24 


8.8 ± 1.4 


16 


32 


8.0 ± 1.2 



Table 1: Dependence of Kss on lattice volume, showing the approach to an infinite volume 
limit. 



where Kss is defined in Eq. (^). Naively, one would convert from lattice to continuum units 
by 

tiss = diffusion constant (lattice units) x {uPLMreY , (20) 
but we have used the thermodynamic corrections derived in the last section, as described in 

The largest 3 lattices are statistically compatible, so we have achieved the large volume 
limit at least at the level of statistics obtained here. The winding number changes on the 
smallest lattice were all immediately followed by a winding number change of opposite sign; 
this might be the system getting almost up to a half integer Ncs state and then turning back, 
and does not represent a permanent change to the underlying vacuum winding number. We 
have only an upper limit for the diffusion constant at that lattice spacing. 

We would like to use this lattice rate to estimate the rate in the physical quantum theory. 
As we have argued, thermodynamic errors, ie errors in ni, are under control and should be 
smaller than our statistical errors. However, the same cannot be said of systematics in 
the dynamics, ie in K2- The problem of relating the classical lattice theory and the real, 
continuum quantum one has recently been studied by Arnold |ll5|. He argues that, in the 



large HTL effect regime, the evolution of infrared magnetic fields is overdamped on time 
scales longer than 1/g'^T, and the strength of the damping sets k,2- The damping occurs 
because E fields set up currents of "hard" modes; the currents propagate; and they enter 
the soft mode equations of motion somewhere else. Hence, the infrared fields at one location 
"feel" the electric fields along the past light cone, due to interactions with hard modes. This 
is the physics of the hard thermal loops. However, the distribution of hard modes on the 
lattice is very anisotropic; so will be the conveying of information by the hard modes; and so 
will be the damping of the infrared fields. Arnold argues that one can make an approximate 
match between classical lattice and continuum quantum values of K2 by taking an angular 
average of the damping strength. There are systematics associated with this, which he 
estimates conservatively as being on order 30%. There is another systematic because the 
lattice system is probably not deep in the strong damping regime; neither is the continuum 
quantum theory at realistic as- 

Using Arnold's proposed match of damping coefficients, we find that Tss for SU(3) with 
six flavors of quarks is (including the factor of 1/2 to go from the diffusion constant to the 
response coefficient) 

r,, = (108 ± 15,tat X 2±i ^y^')alT'' {ml = 2gjT^) . (21) 
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Here I have made a conservative estimate of the systematic errors to be about a factor of 2. 

The ratio Tss/T^g in the formal small a limit has smaller systematics because the problem 
from the anisotropy of the damping coefficient is common to SU(2) and SU(3). Corrections 
to the large HTL limit are probably of different magnitude for the two theories, though 
they are presumably of the same sign. We can estimate them by fitting the SU(2) data at 
different values of /5l to the functional form A/?^^ (leading) + 5/3^^ (correction) and seeing 
how large the extrapolation from = 16 to /S^ = oo is. From the data in [|1^, which 
has the same thermodynamic improvements as here and a topological definition of Ncs, we 
find that the extrapolation is a 21% correction. (The fit is startlingly good: = -42 for 5 
points and 2 fitting parameters.) If the (unknown) difference between the SU(2) and SU(3) 
extrapolations is on order the same size as the SU(2) extrapolation, then taking the ratio of 
/Sl = 16 data and using the size of the extrapolation to give the systematic error gives 

Tss . - - . ^ ^ . / a. \^ mk„f7? 



(10.7±1.53tat±2.3.,.0(f^) 5^- (22) 



ws Kotyj J Tfij^^g 

(The direct ratio of Pl = 16 data is 7.1, but remember that m'jyi^ oc N^.) 

To evaluate this we need and in the dimensionally reduced 3-D theory. Using 
Tc = lOOGeV and as{Mz,MS) = 0.118, we ran to the renormalization point fi = 7.06T 
using the two loop renormalization group equation, and used this value as input in Eq. (146) 
of ||2T|, adapted to six flavor SU(3), to flnd ag^oR- This procedure should minimize two loop 



errors. The result is ag^DR = 0.086. We took from Figure 7 of [^, at Tc = lOOGeV 
and rriH = 70GeV: a^^DR = 1/31.7. Both values and particularly the value for as,DR are 
smaller than we are used to; this is because the dimensional reduction procedure sets the 
coupling roughly to the value at /i(MS) ~ 7T. Using these values, mjj^ = llgl^T^/G, and 



m 



2 _ o^2t^2 
Ds 



^g^T"^, evaluating Eq. (^2]) gives Tss/T^g ~ 1500.| This is much smaller than the 
value we would get using = as{Mz) — 0.118. This large renormalization point dependence 
makes one nervous that subleading contributions in may be non- negligible. 



5 Conclusion 

We have calculated the diffusion constant for SU(3) Chern-Simons number at weak coupling 
(ie, high temperatures) by classical, lattice methods, using a topological deflnition of Ncs, 
and flnd that, at equal values of coupling and hard thermal loops, the diffusion constant 
is larger than in SU(2) by an order of magnitude. Since in the standard model above the 
electroweak phase transition, the strong coupling constant is quite a bit larger than the weak 
coupling constant, the actual ratio of winding number diffusion rates for SU(3) and SU(2) 
is very large. 

Since the hard thermal loop effects induced by the hard lattice modes are different from 
those which would occur for ultrarelativistic particles |]T^, we do not know with certainty 



how to convert the diffusion constant of the classical lattice system into the diffusion constant 
for the physical quantum system, and cannot establish the time constant r with which a 



^Technically we should not just use ag^oR but should also include the correction, Eq. (|l5|). But even 
after taking the fifth power of ag, this correction is less than 10%. This is another example of how well 
behaved the integration over the Aq field turns out to be. 
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chiral quark number is damped by strong sphaleron processes in the plasma. But using the 
estimate, Eq. (^TJ), we get r ~ 80/T. This is slow enough to allow quarks to escape the 
bubble wall before chiral quark number is destroyed, but it is still over 3 orders of magnitude 
faster than the rate at which chiral quark number is converted into baryons through weak 
sphaleron processes. 

We should mention an interesting case where strong sphalerons are less important, which 
is for theories with a stop squark light enough to develop a condensate just before the 
electroweak phase transition. The thermodynamics of this model have been considered 
recently and it apparently provides an especially strong phase transition. Since color is 
broken from SU(3) to SU(2) in the symmetric phase, strong sphalerons only proceed at the 
SU(2) rate for two of the colors and will only erase chiral quark number in the third color 
at an exponentially small rate. Of course, strong processes will mix chiral quark number 
between the two unbroken and one broken color, but the suppression of the rate by a factor of 
10 is significant when the ratio Tss/T^g is important to the final baryon number abundance. 

Acknowledgement s 

I would like to thank Alex Krasnitz, James Hetrick, and Peter Arnold for useful conversa- 
tions or correspondence. This work was supported under NSF contract NSF - PHY96-00258. 

References 

[1] A. Sakharov, JET? Lett. 6 (1967) 24. 

[2] L. McLerran, M. Shaposhnikov, N. Turok, and M. Voloshin, Phys. Lett. B 256 (1991) 
451; A. Cohen, D. Kaplan, and A. Nelson, Phys. Lett. B 263 (1991) 86; M. Dine and 
S. Thomas, Phys. Lett. B 328 (1994) 73. 

[3] A. Nelson, D. Kaplan, and A. Cohen, Nucl. Phys. B 373 (1992) 453; M. Joyce, T. 
Prokopec, and N. Turok, Phys. Rev. D 53 (1996) 2958. 

[4] A. Cohen, D. Kaplan, and A. Nelson, Phys. Lett. B 336 (1994) 41. 

[5] L. McLerran, E. Mottola, and M. Shaposhnikov, Phys. Rev. D 43 (1991) 2027. 

[6] G. Giudice and M. Shaposhnikov, Phys. Lett. B 326 (1994) 118. 

[7] S. Khlebnikov and M. Shaposhnikov, Nucl. Phys. B308 (1988) 885. 

[8] V. Rubakov and M. Shaposhnikov, Phys. Usp. 39 (1996) 461-502, (Usp. Fiz. Nauk 166 
(1996) 493-537.). 

[9] G. D. Moore, Nucl. Phys. B 480 (1996) 657. 

[10] D. Grigorev and V. Rubakov, Nucl. Phys. B 299 (1988) 248. 

[11] J. Ambj0rn, T. Askgaard, H. Porter, and M. Shaposhnikov, Nucl. Phys. B 353 (1991) 
346. 

[12] J. Ambj0rn and A. Krasnitz, Phys. Lett. B 362 (1995) 97. 



14 



[13] G. D. Moore, Nucl. Phys. B 493 (1997) 439. 

[14] P. Arnold, D. Son, and L. Yaffe, Phys. Rev. D 55 (1997) 6264; P. Huet and D. Son, 
Phys. Lett. B393 (1997) 94. 

[15] P. Arnold, Phys. Rev. D 55 (1997) 7781. 

[16] G. D. Moore and N. Turok, PUPT-1681, |hep-ph/ 9703266 . 



[17] D. Bodeker, L. McLerran, and A. Smilga, Phys. Rev. D 52 (1995) 4675. 

[18] James Hetrick, private communication. 

[19] J. Ambj0rn and A. Krasnitz, |hep-ph /9705380| . 

[20] P. Woit, Phys. Rev. Lett. 51 (1983) 638; Nucl. Phys. B 262 (1985) 284. 

[21] K. Kajantie, M. Laine, K. Rummukainen, and M. Shaposhnikov, Nucl. Phys. B 458 
(1996) 90. 

[22] W. Tang and J. Smit, Nucl. Phys. B 482 (1996) 265. 

[23] J. Kogut and L. Susskind, Phys. Rev. D 11 (1975) 395. 

[24] K. Farakos, K. Kajantie, K. Rummukainen, and M. Shaposhnikov, Nucl. Phys. B 442 
(1995) 317. 

[25] J. Mandula and M. Ogilvie, Phys. Lett B 185 (1987) 127. 

[26] G. Parisi and F. Rapuano, Phys. Lett. B 152 (1985) 218. 

[27] D. B5deker, P. John, M. Laine, and M. Schmidt, HD-THEP-96-56, |hep-ph/9612364 . 



15 



